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ABSTRACT 


The ratio of gamma heating per gram of tungsten to gamma heating per gram of water 
was calculated for the case of a thin tungsten detector in a water shield. One-dimensional 
transport theory calculations were used to obtain response polynomials which predict this 
ratio as a function of source energy, shield thickness, and detector thickness. Electron 
transport effects were also examined. Useful results are presented in the form of graphs. 
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SUMMARY 

Two effects must be considered when converting the gamma heating rate measured 
using an aqueous dosimeter to the gamma heating rate in a thin tungsten detector im- 
mersed in a water shield. These are, first, the difference in heating due to the differ- 
ence between the mass energy absorption coefficients of the heavy detector and the aque- 
ous dosimeter, and second, the "skin heating effect. "* The measured heating rates 
must be corrected to account for these effects. Calculated thin detector heating rates 
which are obtained using the buildup factor method also require correction for these ef- 
fects because the buildup factor does not take into account any physical properties of the 
detector. 

One-dimensional photon transport- theory calculations were performed to obtain cor- 
rection factors for these effects for tungsten detectors in a water shield. The results of 
these calculations were used to construct polynomial equations which predict the com- 
bined magnitude of the effects as a function of gamma source energy (from 0. 255 to 
4. 0 MeV), water-shield thickness (from 4. 0 to 44. 0 cm), and tungsten-detector thickness 
(from 0. 00254 to 0. 762 cm). Polynomials were also obtained for the gamma heating 
(W/g) in water and in tungsten so that multienergy sources can be analyzed. 

The behavior of the correction factor (W/g in tungsten divided by W/g in water) 
around the center point (2. 128-MeV source, 24. 0-cm water shield, and 0. 044- cm de- 
tector) was examined using these polynomials. This ratio was also investigated for a 
fission spectrum source. The ratio was less than two for detector thickness greater than 
0. 2 centimeter. The polynomials fit the data with an average absolute error of less than 
10. 2 percent. 

*The "skin heating effect" is the increased gamma heating that occurs near the sur- 
face of the detector. This heating results from the influx of a large quantity of low- 
energy Compton scattered gammas. These gammas are absorbed locally due to the 
large photoelectric cross sections of the heavy element detector. 



INTRODUCTION 


Reactor irradiation experiments may have substantial internal heating by gamma rays 
from the core and surrounding structures. Often it is important to measure the amount of 
gamma heat either for heat-transfer design calculations or to determine effects on tem- 
perature measurements. 

Gamma heating at a given location in water or other light elements can be determined 
using conventional aqueous detectors such as ferrous sulfate dosimeters. However, such 
detectors, although convenient to use, do not give gamma energy spectral information 
which is needed to compute what the heating would be at the same location using a heavy 
element such as tungsten. 

Therefore a conventional way to do design calculations involving heavy elements is 
first to measure the gamma heating using an aqueous dosimeter, obtaining watts per gram 
in water. This watts per gram is multiplied by the mass of an experiment component to 
obtain the heat generation in that component. The result is then multiplied by a correction 
factor for the mass energy absorption coefficient and skin-heating effects. 

There has been a lack of good information on these correction factors. The problem 
has been identified in the literature (refs. 1 to 3), but the choice of the correction factors 
has remained mostly a matter of engineering judgment. It is the purpose of this report to 
show accurate calculated values for these correction factors for the most general case of 
a point source and a shell detector. Also shown are the calculational model and the syn- 
thetic experiment design techniques that were used to select the calculations that were 
performed. The results are presented in the form of polynomial response equations. 
Special cases of interest are shown, and applications of the results are discussed at the 
conclusion of the report. 


SYMBOLS 


A n , A 
o’ n’ nn 

E 

E„ 


F’ 

f 

P 


R 


constant term or coefficient in a response equation 
energy, MeV 

uncoded source energy, MeV 

gamma heating correction factor from response equation 

computed gamma heating correction factor from transport calculations 

fraction of electron energy lost while escaping 

escape probability for electrons 

extrapolated range of electrons in tungsten, cm 
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r radius coordinate in transport model, cm 

S multienergy source, photons/(sec)(MeV) 

T untransformed, uncoded tungsten- detector thickness, cm 
t untransformed, uncoded water-shield thickness, cm 
t* transformed, uncoded water-shield thickness 

Xj coded source energy 

Xg transformed and coded water- shield thickness 
Xg transformed and coded tungsten-detector thickness 
y relative gamma heating, W/g 

9 plane angle corresponding to ft 

\p fraction of electron kinetic energy that is deposited in detector 

ft solid angle 

Subscripts: 

n, nn variable indicators in response equation 

1 number of experiment 

Superscripts: 

W tungsten 

HgO water 


METHOD OF ANALYSIS 

The model selected was a point source of gammas immersed in a spherical water 
shield. The detector consisted of a spherical tungsten shell. A water backstop outside 
of the detector was used with a vacuum boundary condition (see fig. 1). Fluxes were ob- 
tained using DTF-IV, a discrete ordinates computer solution to the Boltzmann equations 
(ref. 4). Scattering cross sections were generated using GAMLEG (ref. 5). Gamma ab- 
sorption cross sections were taken from reference 6. Energy absorption coefficients 
were obtained using data from reference 7 corrected to include energy deposited by 
Compton scattering. To guarantee reliable answers, between 13 and 23 energy groups 
with S|2 Gaussian angular quadrature were used with the Pg scattering approximation. 

The correction factor F was defined as the ratio of the gamma heating in the tung- 
sten detector to the gamma heating in a water detector having the same thickness and 
location (see eqs. (1) and (2)). This ratio was obtained from 15 sets of two companion 
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Figure 1. - Geometrical model used in gamma transport calculations. 


calculations which were the same in all respects except that one calculation of each set 
had a tungsten detector and the other had a water detector. 

To obtain the maximum information possible from a limited amount of calculations, 
the correction factor was studied using multivariate regression analysis. This analysis 
used the calculations as synthetic experiments in an orthogonal composite design (ref. 8). 
This approach produced a response polynomial that may be used to compute the correction 
factor for any values of the independent variables over the range of the variables studied. 
The independent variables selected were source energy, shield thickness (water thickness 
between the source and the tungsten), and detector thickness. A summary of the experi- 
ment statistical design is given in appendix A. 

Because the expected range of correction factors was from 1 to about 30, a least 
squares fit based on minimizing the sum of the squares of the residuals would make the 
fit poor for the smaller values. For this reason the natural logarithm of the correction 
factor was fit instead of fitting the correction factor directly. This fit minimizes the sum 
of squares of the fractional errors at the observed points if the errors are not large (see 
appendix A). 

The orthogonal composite design guarantees that the effects of independent variables 
and their interactions are truly separated (i. e. , their correlation coefficients are zero). 
This gives the most accurate response equation. However, because the levels at which 
two the independent variables must be selected for orthogonality are undesirably spaced, 
new variables were defined which improved the fits. These were obtained by using cubic 
or logarithmic transformations of the originally selected base variables of shield thick- 
ness and detector thickness. These transformations are discussed in appendix A. 

Source energy was not transformed. 


RESULTS 

The correction factor F for monoenergetic source gammas is defined as 
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F(t, T, Ej s 
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y W (t, r, E)r 2 dr dE 
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y 42 (t, r, E)r dr dE 


( 1 ) 


For multienergy sources, the appropriate F factor is defined by 
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F(t, T) = 
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S(E)y W (t, r, E)r 2 dr dE 


100 s>t+T 


If 


H 9 0 9 

S(E)y 6 (t, r, E)r dr dE 


( 2 ) 


where 

S(E) multienergy source, photons/(sec)(MeV) 

y^(t, r, E) gamma heating in tungsten, W/g 
H„0 

y (t, r,E) gamma heating in water, W/g 

The ranges of variables studied were 0.00254- to 0.762-centimeter detector thickness 
(T), 0.255- t 
thickness (t). 

Polynom 

form of these polynomials is 


(T), 0.255- to 4. 0-MeV source energy (E Q ), and 4.0- to 44. 0- centimeter water-shield 

W 

Polynomials were generated for F(t, T, E Q ), y (t, T, E Q ), and y (t, T, E Q ). The 


In K - A Q + AjXj + AgXg + AgXg + A-^X 2 + A^XjXg + AjgXjXg + A 22 x| 


+ A 23 X 2 X 3 +A 33 X 3 


( 3 ) 


where X i5 X„, and X„ are coded and transformed values of the source energy, water- 

W H 2° 

shield thickness, and tungsten- detector thickness, and K is F, y , or y . These 
independent variables are defined by equations (A2) to (A 4) in appendix A. The polyno- 
mial coefficients are shown in table I. 
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TABLE I. - POLYNOMIAL COEFFICIENTS FOR EQUATION (3) 


Coefficient 

K=F(t,T,E 0 ) 

K =y W (t, T, E o ) 

h 9 o 

K = y (t,T,E 0 ) 

A o 

0.8301 

1. 5922 

0. 7640 

A 1 

-.6221 

.2488 

.8697 

A 2 

. 1829 

-.1784 

-.3605 

A 3 

-.8018 

1. 5433 

2.3412 

A 11 

. 5137 

. 0743 

-.4421 

A 12 

-.0084 

.1101 

. 1205 

A 13 

.4169 

. 4173 

-. 0015 

A 22 

1447 

-. 1342 

. 0119 

A 23 

-, 1614 

-.1653 

. 0003 

A 33 

. 0181 

. 0770 

. 0558 



Source energy, E 0 , MeV 


Figure 2. - Gamma heating correction factor F for a 
24-centimeter water shield and a 0.044-centimeter tung- 
sten detector as function of source energy. 


In order to acquire some feel for the behavior of F(t, T, E Q ) when the gamma 
source is monoenergetic, equation (3) was evaluated over the complete range of each in- 
dependent variable around the design center (figs. 2 to 4). These figures may be used 
when the independent variables are approximately those used to draw the curves. How- 
ever, because of the interaction terms, it is advisable to use equation (3) if an accurate 
F factor is desired. The polynomials are valid only over the range of variables studied. 
Hence, extrapolation, especially for thinner detectors, is not recommended. 

W 

The polynomials in y and y can be used to evaluate F for a multienergy 
source (eq. (2)). Equation (2) was approximated for several water- shield and detector 
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Figure 3. - Gamma heating correction factor F for 2. 128-MeV gammas 
and a 0. 044-centimeter tungsten detector as function of shield thick- 
ness. 



Tungsten -detector thickness, T, cm 

Figure 4. - Gamma heating correction factor F for 2. 128-MeV 
gammas and a 24-centimeter water shield as function of 
tungsten-detector thickness. 
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Tungsten -detector thickness, T, cm 

Figured - Gamma heating correction factor F for fission 
spectrum gammas as function of tungsten-detector thick- 
ness T with water-shield thickness t as a parameter. 
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thicknesses by using a five-energy group representation of the Carver prompt plus de- 
layed fission gamma spectrum (ref. 9). These results are shown in figure 5. 


ERROR ANALYSIS 
Uncertainty in the Polynomial Fits 

W 

Table H is a comparison of the calculated F, y , and y values obtained from 
the DTF calculations and the polynomial approximations of these values. The percent 
error at each point and statistical measures of accuracy are also given. The errors at 
the experimental points provide a first-order estimate of the error that can be expected 
for intermediate points calculated using the polynomials (see experiments 16 and 17 shown 
in table H). The statistical measures of fit tend to verify the selection of variables and to 
provide some indication of the overall error and dispersion in the fits. 

If a more exact estimate of the uncertainty is desired, the 95-percent confidence in- 
terval associated with equation (3) may be computed using 

Confidence interval = In F ± 2. 57 a (4) 

(for In F) 


where 

0 = 0. 215s|o. 0667 + 0. 0913(x 2 + X 2 + X 2 ) + 0. 229l(xJ + X* + Xg) 

+ 0. 125o(x 2 X 2 + X 2 x| + X 2 X 2 3 )] 1/2 

and 2. 57 is the two-sided student t for 95-percent confidence with 5 degrees of freedom. 
This means that the true value of In F will lie in an interval whose bounds are ±2. 57 a 
from the polynomial approximation to In F. The constants in this equation were taken 
from the inverted sum-of-cross-products-deviations matrix (ref. 10). Confidence inter- 
vals about F for multienergy sources are complicated by the uncertainty in the integra- 
tion resulting from the energy group scheme selected. This uncertainty was not examined. 
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| 0- 19X8 | 0. 0833 | 0.2126 

l These calculations were performed to test the fit of In F at points not used in determining the polynomial coefficients. 





Use of Numerical Model 


Hopefully some of the inaccuracy due to the numerical model cancels when the results 
of separate solutions are divided to obtain F. There is no general way to assess the ac- 
curacy of an iterative finite difference model. However, because of the large number of 
energy groups, Pg scattering, fine mesh, and tight convergence used, the errors due to 
numerical procedure should be small. 

Errors in Cross Sections 

A detailed discussion of estimated cross section errors is given in reference 6. In 
general, the errors can approach 10 percent below 50 KeV but probably do not exceed 3 
to 5 percent above 100 KeV. 


Electron Transport Effects 

The F values obtained refer to total heat deposition. For the thinnest detectors, 
however, part of the heat is deposited outside of the detector due to transport and escape 
of the electrons. A first-order model shown in appendix B was used to construct a func- 
tion which is defined as the fraction of the electron kinetic energy that is retained in 
the tungsten detector as a function of initial electron energy. In the unperturbed calcula- 
tions (without the tungsten), the net electron transport is negligible. The \f/ function is 
shown in figure 6. 



Figure 6. - Fraction of electron energy retained in tungsten 
detector. 


10 



The correction factor F becomes large only when most of the gamma flux incident 
on the detector is below 1. 0 MeV when the photoelectric effect dominates the gamma heat- 
ing (e. g. , from fig. 3 for a 2. 128-MeV gamma source, F becomes appreciable with a 
few centimeters of water shield thickness). Figure 6 shows that for reactions below 
1 MeV and thicknesses as small as 0.025 centimeter, the correction for electron trans- 
port is small. For thinner detectors the F values should be corrected using the model 
shown in appendix B or a similar treatment. 

USES OF THESE RESULTS 

The polynomials presented here provide a convenient method for computation of heat- 
ing rates in tungsten using measurements made in water. The results have the quoted 
accuracy only for tungsten. However, it is reasonable to assume that other materials 
such as steel, lead, and uranium can be investigated using the same technique. In fact, 
it is suggested that the results given here could be used with reasonable success for other 
combinations of light element shields and heavy element detectors, given that they are 
treated in terms of some "equivalent thicknesses" of water detector and tungsten shield. 

Polynomials such as these could also possibly be used in the place of conventional 
buildup factors in line-of- sight calculations for thin detectors. That is, they provide a 
buildup factor that is a function of the detector in addition to being a function of [it of the 
shield. It should not be difficult to add such a feature to some of the widely used shield- 
ing codes of this type such as QADHD (ref. 11). 

The problem considered here has been examined as a gamma problem without speci- 
fying the type of nuclear reaction that is required to produce the gammas. When these 
results are used to correct gamma heating measurements made in a reactor environment, 
consideration should be given to the fact that secondary gammas from neutron captures in 
the water shield and tungsten detector were not included in this model. Calculation of the 
heating due to secondary gammas requires knowledge of the neutron flux and must be 
treated as a separate problem. When the conversion factors given here are used with 
point kernel integration codes, they can be applied directly. 

As has been noted, the results given here are conservative for very thin detectors 
in that electron losses from these detectors will carry some of the energy away from the 
detector to adjacent structures. The model shown for electron leakage should be useful 
for estimating this effect. Even with corrections for electron transport, however, the 
polynomials should not be used for absolute values of a variable greater than 1. 2154. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 26, 1968, 

120-27-04-54-22. 
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APPENDIX A 


SUMMARY OF EXPERIMENT DESIGN 

To generate the response polynomials, 15 pairs of transport calculations were per- 
formed. Each pair was treated as an ’’experiment. ” The experimental statistical design 
used was an orthogonal composite fractional factorial. Coded levels of independent vari- 
ables for each of the experiments are shown in figure 7. 

Maximum and minimum values were selected for each of the chosen independent vari- 
ables, and the intermediate values were linearly interpolated and coded as shown in refer- 
ence 6. The values used in the DTF calculations are given in table in. 

In order to distribute the calculational information more uniformly over the range of 

-6 -1. -1. 0 0 0 1. 1. 6 Xj level 

0 -1. 1. -6 0 6 -1. 1. 0 X 2 level 



Figure 7. - Experimental design. (6 = ±1.2154 •« the 
value of the highest and lowest levels; numbers 
in blocks are experiment numbers, 1. ) 


TABLE EL - VALUES OF INDEPENDENT 
VARIABLES USED IN DTF 
CALCULATIONS 


Level 

Source energy, 

E o> 

MeV 

Water-shield 

thickness, 

t, a 

cm 

Detector 

thickness, 

T, a 

cm 

-6 

0.255 

4. 00 

0. 00254 

-1 

.587 

12.85 

. 00421 

0 

2.128 

24. 00 

.0440 

1 

3.670 

35.15 

.460 

6 

4.000 

44. 00 

.762 


Require transformation to be orthogonal. 
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water- shield thickness selected for study, the water thicknesses shown in table HI were 
used in the DTF calculations. To maintain the orthogonality of the experiment design, 
the water- shield thickness t was transformed to the new variable t* using 



(Al) 


Also, since the detector thickness effects were anticipated to be largest for the thinnest 
detector, the level spacing dictated by the orthogonal design did not provide information 
in the region of highest importance. For this reason, the log of the detector thickness 
was used as the third variable rather than the dector thickness. 

The equations that define the coded and transformed X values required for evalua- 
tion of the polynomials are 


X- = 

E q - 2.128 

(A2) 

A 1 

1.541 

X 2=( 

t - 24. (A 1 / 3 
v 11.15 / 

(A3) 

111 

CO 

X 

In T + 3. 124 
2.347 

(A4) 


Coded versions of the polynomials have been used here because the resulting response is 
not as sensitive to errors in independent variables as would be the response from polyno- 
mials using uncoded input. 

The least squares fitting procedure minimizes the sum of squares of residuals at ob- 
served points. For the gamma heating correction factor, for example, it is desirable to 
obtain a fit which is as reliable for low values as it is for high ones. To accomplish this, 
the fit was made to In F instead of F. The quantity that was minimized was then 



1=1 



F + dF ~| 2 

F Jl 


(A 5) 


The primes in equation (A 5) denote observed values of F. If the log is expanded for a 
small argument (assuming that dF, the error in F, is small), equation (A 5) becomes 
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15 

2 

~ function minimized by least squares (A6) 

1 

1=1 

This gave a much better fit to the data than the fit of F directly because fractional er- 
rors were being minimized. This reasoning was also used with the polynomials for y W 
HoO 
and y .. 
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APPENDIX B 


DERIVATION OF ijj 

A function was derived which showed, to first order, the fraction of the calculated 
gamma heating that is retained in the tungsten detector. The function is defined as 

/Number of electrons^ /Number of electrons that do escape times the fractionx 

^ _ \ that do not escape / + \ of their energy given up while escaping j_ 

Total number of electrons generated 

= 1 - P(1 - f) (Bl) 

where 

P electron escape probability 

f fraction of electron energy that is lost during escape 

To obtain P, two cases must be considered. First, when R, the extrapolated range of 
the electrons (ref. 12), is greater than or equal to T, the detector thickness, the escape 
probability is given by 


P - 



+ cos fig) dx 


) < 





(B2) 


Figure 8. - Slab model used to obtain P when R> T. 
(The sum of the solid angles, obtained by rotating 
R about the X axis, divided by 47r gives the prob- 
ability of escape for an electron born at X. This is 
integrated over X and divided by T to obtain P 
when R>T.) 
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where the slab model of figure 8 has been used to define 6^ and These angles are 
the plane angles corresponding to the solid angles obtained by rotating the range vector R 
about the X axis. When R < T, a similar integration gives 

P = — (B3) 

2T 

To obtain f, the average distance r through the tungsten travelled by an escaping elec- 
tron is first computed. This is used with the range energy relation (ref. 12) to compute 
the energy required (energy lost) to travel this distance. The fraction of energy lost dur- 
ing escape is this energy divided by the initial energy of the electron. The initial energy 
is assumed to be that of the photon causing the reaction. 

When R > T, 



where use has again been made of the slab model shown in figure 8 and S2 is the solid 
angle corresponding to r = R. 

When R < T, a similar integration gives 

r=- (B5) 

2 

Figure 6 was constructed using equations (B2) to (B5) in equation (Bl). 
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